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ABSTRACT 

We present three-dimensional space velocities of stars selected to be consistent with mem- 
bership in the Virgo stellar substructure. Candidates were selected from SA 103, a single 
40' x 40' field from our proper motion (PM) survey in Kapteyn's Selected Areas (SAs), based 
on the PMs, SDSS photometry, and follow-up spectroscopy of 215 stars. The signature of the 
Virgo substructure is clear in the SDSS color-magnitude diagram (CMD) centered on SA 103, 
and 16 stars are identified that have high Galactocentric-frame radial velocities (Vgsr > 50 km 
s -1 ) and lie near the CMD locus of Virgo. The implied distance to the Virgo substructure from 
the candidates is 14 ±3 kpc. We derive mean kinematics from these 16 stars, finding a radial ve- 
locity Vgsr = 153 ±22 km s -1 and proper motions (fi a cos5,/j,$) = (-5.24, -0.91) ±(0.43, 0.46) 
mas yr -1 . From the mean kinematics of these members, we determine that the Virgo progenitor 
was on an eccentric (e ~ 0.8) orbit that recently passed near the Galactic center (pericentric 
distance R p ~ 6 kpc). This destructive orbit is consistent with the idea that the substructure(s) 
in Virgo originated in the tidal disruption of a Milky Way satellite. Af-body simulations sug- 
gest that the entire cloud-like Virgo substructure (encompassing the "Virgo Overdensity" and 
the "Virgo Stellar Stream") is likely the tidal debris remnant from a recently-disrupted massive 
(~ 10 9 M Q ) dwarf galaxy. The model also suggests that some other known stellar overdensities 
in the Milky Way halo (e.g., the Pisces Overdensity and debris near NGC 2419 and SEGUE 1) 
are explained by the disruption of the Virgo progenitor. 

Subject headings: Galaxy structure - Galaxy: kinematics and dynamics - Galaxy: stellar con- 
tent - Proper motions - stars: kinematics and dynamics - galaxies: dwarf 
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1. Introduction 



Galaxy formation is currently thought to proceed via the hierarchical merging of numerous smaller 
subst ructures, which are assimilated into the host galaxies by tidal disruption and , ultimately, dissolved 
(e.g.. ljohnstonlll998l : IXbadi et allbooa: bullock & Johnston] bood : Font et allboOQ) . This is by no means 
a new idea, as it was suggested by ISearle & Zinnl (119781) that late infall after the collapse of the initial 
protogalaxy w as responsible for the apparent age spr ead among Milky Way (MW) globular clusters, and 
later (by, e.g., lMajewskilll993l : LMajewski et all 1 19961 ) that the build-up of the Galactic halo proceeds via 
agglomeration of dwarf galaxies. We have now seen this process in a ction via the nu merous stellar tidal 
streams that have been discovered coursing through the MW's halo ( see iGrillmainBO 1 Ol for a summary) and 
various other s tellar substructures in the form of "clouds" of stars (the properties of known structures are 
summarized in iRocha-Pintol I2010I) . In the era of large-area, deep photometric and spectroscopic surveys, 
we are now entering the regime where not only narr ow, kinematicaily-col d stellar streams can be found, but 
also the clouds and shells of debris seen in, e.g., the ljohnston et al. ( 2008 ) models of hierarchical formation 
of Milky Way-like galaxies. The menagerie of morphological structures predicted by these models are now 



being seen with deep imaging around external galaxies (e.g.. IChonis et al.l 1201 ll ; iMartfnez-Delgado et al. 



2010). Because of our perspective within the Galaxy, such features are more difficult to see in context 
and characte r ize. H owever, a look at the variety of morphological types of known structures discussed in 
Rocha-PintoN20lrt) shows that large photometric surveys have begun finding the predicted remnants of halo 
assembly. The kinematical characterization of the clouds and shells that is needed to assess their origins 
and discern the nature of their progenitors has been more elusive, since identification of stellar members of 
such low surface-brightness features is difficult, and derivation of proper motions for such distant structures 
is challenging. The particular one of these apparent cloud-like stellar structures upon which we focus this 
work is the overdensity (or, perhaps, collection of stellar excesses) in the constellation Virgo. 

The first recogni tion of a stellar overdensity in the halo near the Virgo constellation was the finding 
by Ivivas et all (feoOll ) of a cluster of 5 RR Lyrae stars within 195° < RA < 200°, at a dist a nce o f ~20 
kpc in the first 100 deg 2 results from the QUEST survey. Shortly thereafter, INewberg et all (|2002r ) used 
F-type turnoff starcounts on the celestial equator to identify stellar overdensities, including one at 180° < 
RA < 195°, Dec = 0° that was dubbed S297+63-20.0 (the naming scheme denotes that the feature is at 
(l,b)~ (297° , +63 ° ) with turnoff magnitude g ~ 20). This feature, which is near the same region of sky as the 
Vivas et al.1 (120011) RR Lyrae, shows a clear main sequence extending downward from g ~ 19, corresponding 
to a distan ce (assuming an ol d, metal-poor population) of ~ 20 kpc. Since the distance is also similar to 
that of the 



Vivas et al. 



(2001) stars, this large-area overdensity might be associated with the RR Lyrae. 
Newberg et al.l (|2002l) also noted that the color of the turnoff stars in S297+63-20.0 is intermediate between 
that of Sagittarius tidal debris and the typical MW halo turnoff, suggesting that the Virgo feature is related 
to neither of those structures. 

Further result s from the QUEST survey, which studied a 2.3-degree wide strip centered at Dec ~ -1°, 
were presented in lVivas & Zinnl (12003b . In the 380 deg 2 portion of the QUEST survey discussed, the authors 
find 497 RR Lyrae between 4-50 kpc. Many of the distant RR Lyrae are Sgr debris, but the authors identify 
a distinct clump at mean magnitude of (Vq) = 16.9 ±0.2, corresponding to a distance of ~ 19 kpc. Because 
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of its location at (RA) ~ 186°, this grouping was dubbed the "12.4-hour clump". Based on the RR Lyrae 
periods and period-amplitude distributions, it appears that the clump contains RR Lyrae of both Oosterhoff 
classes, making them more typical of dwarf galaxy populations than globular clusters. This 12.4-hour 
clump is within the same region and at the s ame distance as S297+ 63-20.0, which suggests an association. 
A "fluff" of overdensity was also seen in the lMajewski et all (120031) study of 2MASS-selected M-giant stars 
at the same position (150° < RA < 210°) as S297+63-20.0 and the 12.4-hour clump, though it was suggested 
by those authors to be part of the "descending, foreshortened loop" of the Sagittarius leading tidal tail. This 
idea was explored further by IMartmez-Delgado et all (120041) . who develop ed a model of the Sgr disruption 
that pr oduced debris a t roug hly the distance, position, and velocity of the IVivas & Zinnl (12003b RR Lyrae 
and the iNewberg et al.l (120021) main sequence detections. 



.T 



The first spectroscopic study of any of these overdensities came when Duffau et al. (2006) observed 18 
RR Lyrae from the QUEST survey. Their study found that a clump of stars centered at a ~ 186° (±5°) and 
^Sun = 19.6 kpc is also clustered in velocity. The six stars making up this clump have a mean (Vgsr) = 99.8 
km s , and a dispersion of 17.3 km s -1 . This clump of stars was dubbed the "Virgo Stellar Stream (VSS)". 
Duffau et al. then expanded the sample to a larger area and RVs within ~ 3a of the peak, and identified 
13 stars with (Vgsr) = 85 km s" 1 with ~ 25-30 km s -1 dispersion. These stars are spread over a distance 
range from /?s U n = 16-20 kpc (with one at ~ 24 kpc); thus, the VSS must be large and diffuse (or contain 
multiple overdensities). The metallicity of the 7 RRab stars selected to be members is ([Fe/H]) = -1.86, 
with <7[Fe/H] = 0-4. The large spread in conjunction with the lack of an obvious intermediate-age population 
in the SDSS CMDs is taken as evidence of a low-luminosity, old, metal-poor dwarf galaxy progenitor for 
the VSS. 

Furt her detections of halo overdensities from the QUEST RR Lyrae survey's first catalog were pre- 
sented by IVivas & Zinnl (120061) . who found that after removing the contribution of their best-fitting smooth 
halo model, the Virgo structure is visible as a peak at RA~ 186° at a distance of ~ 17 kpc (spanning 12-20 
kpc betw een 175° <RA< 2 05° ). Velocity substructures were seen among bright (V < 16.1) QUEST RR 
Lyrae by IVivas et al.l (12008b . including a grouping at 10.5 < D < 12 kpc with (Vgsr) = 215 km s -1 and 25 
km s -1 dispersion (plus some stars at slightly lower velocities that may also be associated with the struc- 
ture). These stars extend the Virgo substructure to more nearby distances than previously known. The mean 
metallicity of six of these bright RR Lyrae is ([Fe/ H]} = —1.55 , with opwH] = 0.15 dex. Further spectro- 
scopic data of QUEST RR Lyrae was reported by Duffau et al.l (|2010b . This study found a VSS grouping 
with 23 members, having a mean RV of (Vgsr) = 129 km s _1 . There are additional stars in a group at higher 
(200 < Vgsr < 250 km s" 1 ) velocity, bu t slightly more nearby (9-15 kpc). These, together with the 10-12 
kpc ~ 200 km s" 1 stars from I Vivas et al.l (120081) suggest a velocity gradient with distance in the VSS region. 



Table 1 : Previous detections of stellar overdensities in Virgo. 



source 


RA 


Dec 


distance 


Vgsr 


[Fe/H] 


notes 




(degrees) 


(degrees) 


(kpc) 


(km s" 1 ) 


(dex) 




Vivas et al. (2001) 


[195,200] 


[-2.3,0] 


~20 






5 QUEST RR Lyrae 


Newberg et al. (2002) 


190.2 


0.3 


~20 






SDSS F-turnoff stars 


Vivas & Zinn (2003) 


186 


-1 


19 






QUEST RR Lyrae (" 12.4-hour clump") 


Duffau et al. (2006) 


[175,205] 


[-2.5,0.0] 


[16,20] 


99.8 


-1.86 


spectra of QUEST RR Lyrae ("Virgo Stellar Stream") 


Vivas & Zinn (2006) 


189 


-0.8 


17 ([12,20]) 






QUEST RR Lyrae 


Newberg et al. (2007) 


191.2, 185.6 


-7.8, -0.1 


18 


130±10 




SDSS F-turnoff stars 


Juric et al. (2008) 


[160,210] 


[-5,15] 


6-20 




<-i.o 


SDSS "tomography" ("Virgo Overdensity") 


Keller et al. (2008) 


[180,230] 


[-12,-2] 


[16,19] 






SEKBO RR Lyrae 


Vivas et al. (2008) 


[170,210] 


[-2.3,0.0] 


[8,12] 


215 


-1.55 


spectra of bright (V < 16.1) QUEST RR Lyrae 


Keller et al. (2009) 


[180,210] 


[-15,0] 


[17,20] 






SEKBO RR Lyrae with SDSS photometry 


Prior et al. (2009) 


[170,210] 


[-14,0] 


[16,22] 


127 ± 10 a 


-1.95 ±0.1 


spectra of SEKBO RR Lyrae 1 


Brink et al. (2010) 


185.9 


-1.0 




137 b 




spectroscopy WFC on INT 2.5-meter 


Duffau et al. (2010) 


[180,210] 


[-4,0] 


[7,21] 


129 




"VSS" group of QUEST RR Lyrae 


Duffau et al. (2010) 


[180,210] 


[-4,0] 


[7,18] 


~220 




higher-velocity group of QUEST RR Lyrae 


Keller (2010) 


[180,210] 


[-2,6] 


[12,18] 






SDSS subgiants 


Casey et al. (2012) 


[180,194] 


[-3,0] 


~20 


~ 130 


-2.0 


spectra; Feature B 


Casey et al. (2012) 


[180,194] 


[-3,0] 


~20 


~220 


-1.2 


spectra; Feature C 



Values in brackets represent ranges in which those particular studies found Virgo substructure members. 

"Note that this study contained some stars at high (Vgsr > 180 km s -1 ) velocity that were not identified by the authors as VSS 
members. 

*We note that the authors' estimate of V G sr = 137 km s -1 uses only stars with Vgsr < 170 km s -1 ; using all stars with 100 < V G sr < 350 
km s" 1 yields a median Vqsr = 152 km s -1 . 
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The overdensity of R R Lyme in Virgo has also been studied by the SEKBO survey, which spans ~ 10° 
about the ecliptic. iKeller et al.1 ( 120081) show that the Virgo overdensity is visible among SEKBO RR Lyrae 
at R ~ 20 kpc be tween 180° <RA< 230°. Photometric and spectroscopic follow-up of SEKBO stars by 



Prior et al.l (120091) uncovered a handful of additional VSS stars between 16-22 kpc at mean Vgsr = 127 ± 10 
km s -1 and dispersion of ~ 27 km s -1 . These stars have ([Fe/H]) = -1.95 ±0.1, with scatter of 0.4 dex 
(note that an additional star at Vgsr = 193 km s -1 and [Fe/H] = -1.34 was not included as a member in their 
sample). [Prior et al.l (|2009i) also calculate a number density of excess stars and use this to find that in a 
~ 760 deg 2 area centered on (a, 5) = (186°, -4°) the total luminosity is My = -11.9 (or My = -10.1 using 
values more closely resembling the IJuric et al.l 120081 distances). The authors claim that becuase VSS stars 
have been found between 12-24 kpc with a large metallicity spread, we must be looking at the remains of a 
disrupted dwarf spheroidal galaxy. We also note that this study extended the VSS to much lower latitudes 
of Z?~45° (<5~-15°). 

An extens i ve exa mination of the Virgo substructure was enabled by SDSS Data Release 5 (DR5) data. 
Newberg et al.l (|2007l) studied Virgo using SDSS starcounts of BHB stars, blue stragglers, and F-turnoff 
stars as well as SEGUE radial velocities of F-type turnoff stars in the region of the Virgo overdensity. The 
authors showed using BHB stars and blue stragglers within 15 kpc of the Sgr orbital plane tha t Sgr debris in 
this region does not pass through the position of S297+63-20.5 identified by iNewberg et al.l (12002b . ruling 
out the association of the Virgo structure with leading-arm Sgr debris. Since the SDSS footprint by DR5 
covered the entire north Galactic cap (NGC), Newberg et al. were able to make starcount maps of the NGC 
using photometrically-selected F-turnoff stars. These maps showed that the Virg o (S297 +63-20.5) struc ture 
is prominent mostly at magnitudes 20 < g < 21, spanning at least ~ 15° across. INewberg et al.l (|2007|) use 
SEGUE radial velocities of F-turnoff stars in two plates at (l,b) = (300°, 55°) and (288°, 62°) to search for 
a Virgo RV signature. A clear excess of F stars is seen at (Vgsr) = 130 ± 10 km s -1 with stars in the 30 km 
s -1 wide bins adjacent to the peak, suggesting a large velocity dispersion within the structure. 



Brink et al.l ((2010) undertook a spectroscopic survey of 111 faint (18.0 < go < 21.8), blue (-0.1 < 



(g~ r )o < 0.65) turnoff candidates centered on the VSS, the majority of which seem to be Sagittarius leading 
arm stars. However, there is a b road peak at (Vgsr) = 137 km s -1 with dispersion 22 km s -1 that is consistent 



with the INewberg et 



al 



(2007) V SS/VO D detections from F-turnoff stars. We note also that the mean and 
dispersion quoted by iBrink et al.l (|2010Q was measured using only stars between 105 < Vgsr < 170 km s , 



but there is a long tail of stars at higher velocities in their sample. 



The final spectroscopic study we mention is that of ICasey et al.l (120121) . These authors observed a 
sample of 178 K-giant stars in four separate fields between 180 < RA < 195° and -3 < Dec < 0°. In 
addition to a clear signal from the Sagittarius tidal tail at Vgsr < km s" 1 , they found a significant excess 
of high- velocity (Vgsr ^ 70 km s -1 ) stars. The authors suggest that these ~ 20-kpc-distant stars may be two 
separate features (denoted "Feature B" and "Feature C") with 70 < Vgsr ^5 250 km s . 



Another SDSS study by lJuric et al.l d2008l) explored the photometric properties of the overdensity in 
Virgo. The prominent factor of ~ 2 stellar excess spanning over 1000 deg 2 in this study was dubbed the 
"Virgo overdensity (VOD)" by the authors. The VOD extends from the lower limit of the SDSS footprint 
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at Z ~ 6 kpc above the Galactic plane to the outer limit of their study at Z ~ 20 kpc, suggesting that it 
may extend even farther. The feature appears to narrow at higher Z, which suggests that it may extend 
further into (or beyond) the plane. Assuming a distance of 10 kpc and 1000 deg 2 area, IJuric et all (|2008) 
estimated a lower limit on the total luminosity ofL,=0.1 x 1O 6 L , orM r = -7.8. The color-color distribution 
of VO D stars suggests that t he VOD is more metal-poor than the thick disk. However, Juric et al. (as 
well as iMajewski et al.1 120031 previously) see the VOD in 2MASS -selected M giants, so there must also 
be a relatively metal-rich population associated with the VOD as well. Again, because of its low surface 
brightness, large extent, well-defined outline, and low metallicity, the authors argue that the VOD is the 
result of a low-metallicity dwarf galaxy merging with the Milky Way. 



An et al. 



(2009) explored 



Additional work using SDSS-derived photometric metallicity estimates by 
the entire NGC region. The VOD in the region between 270° < / < 330° and 60° < b < 70° was found by 
these authors to have median metallicity of [Fe/H]=-2.0±0.1. Note, however, that this is for all dwarf stars 
in this region, only some of which are VOD members. An et al. claim that the fact that the median [Fe/H] 
derived is the same as the halo metallicity in the symmetric area of the sky means that the additional VOD 
stars are halo-like. 



Using SDSS photometrically-selected sub-giants, iKellerl (120 lOi) mapped the north Galactic cap, an d 
resolved the VOD into a number of discrete overdensities of smaller area than that of Uuric et al.1 ((2008). 
Keller et all ((2009) found a similar result using combined SEKBO RR Lyrae data and SDSS photometry. 
The VOD appears to break up into at least three distinc t clumps in that st udy. This may suggest significant 
inhomogeneity across the larger ~ 1000 deg 2 feature of ljuric et all (12008b . 



A summary of the numerous detections of stellar overdensities in Virgo discussed in this section is 
given in Table [T] From this table, it is obvious that there is a large complex of structures in this part of the 
sky detected using a variety of stellar tracers, that (based on their similar distances and velocities) may be 
related to each other. Each of the surveys presents a data set limited either in its spatial coverage or based 
on a modest-sized sample of spectra, making interpretation difficult. 

Understanding the nature of the Virgo overdensitiesQ and their possible inter-relationship will ultimately 
require measurement of the kinematics and/or chemistry of large numbers of stars in the large area where 
overdensities have been detected. In contrast, we focus on a small region of sky, but where we can take 
advantage of precise measure ments of three-dimensi onal kinematics of selected Virgo candidates. This 
study is a follow-up to that of ICasetti-Dinescu et al.1 (2009; hereafter CD09), which presented the orbit 
of a single RR Lyrae consistent with membership in a Virgo substructure. This star had radial velocity 
measurements from SDSS and proper motion measureme nts from photographic plate data derived from the 
survey of Kapteyn's Selected Areas described initially by lCasetti-Dinescu et al.1 (|2006l ). The 3-D kinematics 
of this particular star place it (and, by inference, the substructure of which it is part) on a highly eccentric 
orbit, having just passed its pericenter. The authors considered possible positional associations with Milky 



'Because the relationship between the various structures including the "Virgo Overdensity (VOD)" and "Virgo Stellar Stream 
(VSS)" in the Virgo region is unclear, we will simply refer to "the Virgo stellar (sub)structure" throughout this work. We will show 
later in this work that all of these structures are likely related, making this semantic approach reasonable. 
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Way satellites beyond Galactic radii of 50 kpc, and found that the unusual globular cluster NGC 2419 is 
within 13 kpc of the derived orbit. 

The present study revisits the result of CD09 using spectra for additional stars from the same proper- 
motion catalog used in the CD09 study. We spectroscopically identify 16 Virgo members, from which 
we derive mean proper motions and RVs. From these, we derive an orbit for the Virgo substructure, and 
use this orbit to produce Af-body models based on the measured kinematics. The paper is structured as 
follows: In Section 2 we discuss the data used for this study, which include proper motions from the Selected 
Areas survey, SDSS photometry, and low-resolution spectroscopy obtained on three separate observing runs. 
Section 3 outlines the methods we used to select Virgo substructure members, and compares our data sets 
to SDSS data from the surrounding regions of sky. The criteria adopted to identify Virgo structure members 
were based on radial velocity, color-magnitude location, proper motion, and to some extent metallicity, 
yielding a final sample of 16 Virgo members. We follow this, in Section 4, with analysis of the kinematics 
derived from the identified members. We derive an orbit based on the mean motion we measure, and use 
this to generate an Af-body simulation that roughly reproduces the (still poorly-known) characteristics of 
the Virgo substructure(s). In Section 5, we search for positional and kinematical associations of known 
Milky Way satellites and substructures with our derived orbit for Virgo debris, and discuss a few possible 
associations. In particular, our data and subsequen t simulations suggest plausible associations of the Virgo 
debris with NGC 2419, the Pisces Overdensity, the Newberg et al.l (12003!) 90-kpc "Sgr debris", and a stellar 
overdensity in the field of view of the ultrafaint dwarf galaxy SEGUE 1 . We conclude by arguing that we 
have convincingly shown that all or most of the known Virgo overdensities are plausibly part of a single 
disrupted dwarf galaxy, and that a number of other Milky Way substructures may be associated with this 
same progenitor. 



2. The Data 



2.1. Photometry and Proper Motions 



Proper motions used in this study come from the survey described in Caset ti-Dinescu et all ([2006, 
hereafter CD06); we provide a brief introduction here, and refer the reader to ICD06l for details. The survey 
consists of ~50 fields of 40' x 40 / sp aced ~ 15 degrees apart in strips along the celestial equator (i.e. ,6 = 0°) 
and 5 = ±15° (see Fig. 1 in |CD06I for a map). These fields are the "Selected Areas (SAs)" designated by 



Kapteynl ( 1906) for his proposed comprehensive study of the Milky Way. The particular field of view in 
this study, SA 103, was also the focus of a previous publication (|CD09|) . SA 103 is centered at (ck, 5)aooo = 
(178.765°, -0.562°), or (l,b) = (274. 56°, 59. 17°). The photographic plates making up our earliest epoch 
for proper motions were ori ginally taken between 1909-1912 (analysis of photometry from these plates 
appears in ISeares et al.Nl93O0 as part of the Kapteyn survey at the Mt. Wilson 60-inch. Matching plates (in 



field of view, depth, and plate scale) were taken by S. Majewski between 1996-1998 at the Las Campanas 
du Pont 2.5-meter telescope to provide a ~ 90-year baseline for measurement of proper motions. The 
Mt. Wilson 60-inch plates have plate scale 27" 12 mm -1 , and the plate scale for du Pont 2.5-meter data is 
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10"92 mm -1 , with both having a field of view of 40' x 40'. All of these plates were digitized using the 
Yale PDS microdensitometer. The data were supplemented by intermediate-epoch plates from the Palomar 
Observatory Sky Survey (POSS-I), taken in the early 1950s, to extend the limiting magnitude of our proper 
motions beyond the depth provided by the Mt. Wilson plates. Proper motions were measured by comparing 
the positions of objects on 1 Mt. Wilson plate, 4 POSS-I plates, and 2 du Pont plates. Each of the Mt. Wilson 
and du Pont plates contain a long and short exposure offset spatially, providing 2 measurements per plate. 
Two separate scans of the POSS-I plates were analyzed; one from the PMM measuring engine at USNO, 
and the other from the PDS microdensitometer at STScI. Thus there is a total of 14 possible measurements 
from the 7 plates in SA 103. The correction from relative t o absol ute proper motions was done using four 



QSOs and 302 background galaxies in the field of view (see CD09 for details), from which we derived zero 
points of fi a cos 5 = 3.89 ±0.27 mas yr -1 andag = 2.31 ±0.30 mas yr ; these were applied to the data to 
place the proper motions in an absolute frameo QSOs are taken from SDSS QSO catalogs, and galaxies are 
identified in the plate scans. Proper motion uncertainties range from ~ 1 mas yr -1 per star at the bright end 
(i.e., g < 18) to 2-3 mas yr" 1 at the faint (g ~ 20-21), limiting magnitude of the plates. 

Rather than relying on photographic photometry, we took advantage of the fact that high-quality pho- 
tometry from SDSS DR7 was available for the SA 103 are a. The photometric d ata we have used throughout 



this work are corrected for reddening and extinction (from ISchlegel et al.ll 19981. the mean reddening is E(B- 
V) = 0.026 along this line of sight). By comparison of photometric catalogs to the much deeper SDSS data, 
we have found that the proper motion catalogs are 86% complete for point sources at g = 20.0, falling to 
15% completeness at g = 21.0. The limiting magnitude set by the shallow Mt. Wilson plates (i.e., for stars 
with ~ 90 yr baseline) is g ~ 20.0, with 90% completeness at g = 17.0. A color magnitude diagram is given 
in the upper left panel of Figured] for all stars with measured proper motions in SA 103 . For comparison, 
the middle panel shows the predictions of the Besancon Galaxy model ( Robin et al 1 l2003h for the same field 



of view, and the right panel shows all stars from SDSS DR8 within the same field of view. In the DR8 panel, 
an upper main sequence is clearly visible, with a turnoff at g ~ 19.0, g -r ~ 0.3 - this likely corresponds 



to the MSTO of the Virgo structure identified by INewberg et al.l (|2002). The same feature is visible in our 



SA 103 data, and is not represented in the Besancon model, which bolsters our assertion that this turnoff is 
due to a substructure in the Milky Way halo. To further confirm that this turnoff structure is not expected 
among Galactic stellar populations in the SA 103 field of view, we compare SDSS data in SA 103 to SDSS 
CMDs at 90° intervals in longitude and approximately the same latitude (b = 60°). Assuming the Milky Way 
has an axisymmetric stellar halo, these four high-latitude fields seen in Figure [Qshouid exhibit the same fea- 
tures among faint, blue halo stars. Indeed, the three fields at / = 0°,90°, and 180° (the bottom panels in 
Figure QJ look similar to each other, but the SA 103 field of view (the upper right panel) shows a clear main 
sequence extending downward from g ~ 19, g-r ~ 0.3. Within the overplotted box at 0.15 < g-r < 0.5 and 
18.5 < g < 20.5, there are a total of 93, 130, 85, and 183 stars for / = 0,90, 180,274° respectively. These 
correspond to fractions of 0.123, 0.125, 0.106, and 0.142 of the total number of stars above the dashed line 
(i.e., with g < 0.7(g-r) + 20.8) in each panel of Figure Q] We thus conclude that the apparent clumping of 



2 Note that these differ from the zero points in 



CD09, but are consistent within the errors, and thus make little difference in our 



final conclusions. The difference between the measurements is in the sample of QSOs used in the two, separate measurements. 
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stars at g > 18.5 and 0.15 < g—r < 0.5 in our SA 103 proper motion data near / = 274° is due to a significant 
stellar excess that is present along this line of sight. Based on these starcounts we surmise that > 30% of 
the stars in this CMD-selection region are part of the Virgo substructure. 

2.2. Radial Velocities 

In Figure Q] there is no obvious red giant branch associated with the MSTO feature, meaning that 
spectroscopic follow-up aimed at efficiently identifying and measuring the kinematics of Virgo substructure 
members should ideally focus on the high-density turnoff region. However, many of the spectra we observed 
for this study were taken as backup targets during observing time scheduled for other projects. Thus the 
spectroscopic targets, shown as large colored circles in Figure |2] are not optimally selected and tend to 
sample brighter magnitudes. The box overplotted in Figure |2](0. 15 < g—r < 0.5, 18.5 < g < 20.5; the same 
as in FigureQ} surrounds the apparent MSTO of the Virgo structure, and thus contains the best spectroscopic 
candidates. Of the total of 215 stars for which we obtained spectra, 53 are within this box, with the majority 
of the remaining spectra at brighter magnitudes. Many of the stars brighter than g ~ 18.5 were selected in 
the region immediately surrounding the RR Lyrae star (g = 17. 05, g-r = 0.24) on which the CD09 study 
focused, with the hope of identifying additional AGB/HB stars in the Virgo structure. The remainder of 
bright targets are redward of the thin/thick disk MSTO, where Virgo RGB candidates would be expected to 
reside. 

Spectra of candidates were obtained on three separate observing runs. The first of these, in June 2009, 
used the Hydra multifiber spectrograph on the WIYN 3.5-meter telescope]! A single setup consisting of 70 
targets was observed at R ~ 1500 (3.35 A at A = 5200 A ) using the 600-line ("600@ 10.1") grating centered 
at A ~ 5700 A (providing wavelength coverage of 4300-7100 A). The remaining candidates were observed 
as backup targets during runs at Cerro Tololo Inter-American Observatory in February 2010 and January 
2011 using the Hydra multifiber spectrograph on the 4-meter Blanco telescope. The observing time for 
these runs was granted by NOAO to search for QSOs to use as astrometric reference sources in the field of 
view of the Carina dwarf spheroidal - results from that study will be presented in a future work. We observed 
at low resolution (R ~ 850), using the KPGL2 grating in first order centered at A ~ 6000 A. This setup was 
chosen to provide spectra that cover a broad (A « 3400- 8200 A) wavelength range for finding the emission 
lines for QSOs within a large redshift window. Such spectra are not ideal for stellar kinematics, providing 
radial velocities of ~ 10-20 km s -1 accuracy (compared to 5-10 km s -1 with the WIYN+Hydra setup; both 
estimates assume cross-correlations accurate to ~ ^ of a resolution element). However, in both the 2010 
and 201 1 runs, our observations were scheduled in a period when the Carina dSph was unobservable at the 
end of the night, so we observed SA 103 to fill extra time. In February 2010, we observed 3 configurations 
totaling 188 stars that had spectra of sufficient signal-to-noise (S/N > 10) to derive velocities. In January 
2011 we observed two additional SA 103 configurations with the same spectrograph setup, totaling 162 



3 The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, Yale University, and the 
National Optical Astronomy Observatory. 
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Fig. 1. — Color magnitude diagrams comparing all stars with measured proper motions in our SA 103 
data (upper left panel) at (l,b) = (274.6°, 59.2°) with output of the Besancon model along the same line 
of sight (upper middle panel) and SDSS CMDs from one-degree fields at b ~ 60° at ~ 90° intervals in 
Galactic longitude (the upper-right and lower three panels, with l,b as labeled; the upper right field is SDSS 
data within the SA 103 field of view). For SA 103 we show each individual star, while all of the other 
panels are binned Hess diagrams, with the grayscale density values given by the colorbar in each panel. 
The dashed line in each panel, defined by g < 0.7(g-r) + 20.8, roughly represents the magnitude limit 
of the SA 103 data, and is given for reference. Assuming the Galactic halo is spherically symmetric at 
such high latitudes, the starcounts of halo stars within the box overplotted on each panel should be roughly 
equal at each longitude. We find a total of 183 stars within the CMD-selection box in the upper-right panel 
(corresponding to SA 103), well above the totals of 93, 130, and 85 that are found in each of the fields at 
/ = 0,90, 180°, respectively. When we instead consider the fraction of the total number of stars above the 
dashed line that are within the boxed region, we find 0.123, 0.125, 0.106, and 0.142 for / = 0,90, 180, and 
274° , respectively. We posit that this excess over the number of stars in symmetric fields (all selected from 
data of equal quality and depth) signals the presence of a significant overdensity (comprising ~ 30% of the 
stars in the boxed region) in SA 103. 
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Fig. 2. — Color-magnitude diagram of all stars with measured proper motions in SA 103 shown as small 
black dots. Stars with spectroscopic data are highlighted by larger circles color-coded by Vqsr- Brighter 
colors correspond to higher Vgsr- The boxed region corresponds to the same Virgo turnoff selection area as 
in Figured] Note that a larger fraction of the light-colored (i.e., high- velocity) stars is present in the area of 
the Virgo turnoff than outside this box. 

stars with usable spectra. Many stars were targeted on multiple configurations to allow cross-checking of 
the radial velocities; thus the total number of stars with measured velocities in SA 103 is 215. 

For all observing, the target configurations were observed multiple times to enable cosmic -ray removal. 
The initial pre-processing of 2-D spectra used the CCDRED package in IRAF0 One-dimensional spectra 
were extracted from the summed 2-D frames using the DOHYDRA routines (also in IRAF) for multifiber 
data reduction. Dispersion solutions were determined using 30-35 emission lines from arc lamp exposures 



4 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation. 



-12- 



taken at each Hydra configuration. Bright radial velocity standards of various spectral types from F through 
K (both dwarfs and giants) were observed on each run. Each of these was observed through multiple fibers, 
to yield numerous cross-correlation spectra for RV measurement. Cross-correlation of these standard spectra 
against each other was done using the IRAF FXCOR package; we found that measured velocities of the RV 
standards typically agreed with published values to within 1-2 km s -1 for the June 2009 WIYN+Hydra 
run, and to < 4 km s -1 for the lower-resolution Feb. 2010 and Jan. 2011 CTIO+Hydra observing. The 
radial velocities for program stars were then derived by cross-correlating the object spectra against all of the 
standard stars observed on the same run on which the target star was observed. 



Radia l velocity uncertain ties were derived using the lVogt et all (11995b method, as desc ribed in lMunoz et al 
d2006h and barlin etaD d2012h . Because the Tonry-Davis ratio (TDR; lTonry & Davislll979h scales with S/N, 
we were able to use the dependence of the TDR on S/N mapped via multiple observations of each standard 
star (which of course have varying S/N) to derive individual RV uncertainties from the measured TDR of 
each spectrum. Typical RV uncertainties for individual measurements are oy « 10-20 km s -1 for spectra 
having S/N >15-20 per Angstrom. For the 100 stars (out of 215 total) that were observed multiple times, 
the final radial velocities represent the error- weighted mean values of individual measurements. We note 
that the RV uncertainties based on the TDR method may be underestimated; comparisons of the veloci- 
ties of stars observed multiple times show standard deviations of 25-30 km s . As this includes all RV 
measurements (including those from low S/N spectra), we estimate that typical uncertainties are ~ 20 km 



2.3. Metallicities 



To estimate metallicities from the l ow-resolution spectra used in th is study, we relied on the well- 
calibrated Lick spectral indices (see, e.g.. lWorthey et al.lll994l ; lFriellll987r) . Details of the software pipeline 
produced by our group for this puipose will appear in a forthcoming paper (Carlin et al. 2012, in prep.). 
The program derives [Fe/H] bas ed on fits of the beha vior of [Fe/H] as a function of eight Lick Fe indices 
and the H/3 index as defined by Worthe v et all (1 1994) . The fits are calib rated using the known [Fe/H] and 



spectral indices in the atlas of Schiavon ( 2007 . based on the spectra of Jones 1998|). This code has been 
shown to yield [Fe/H] measurements of ~ 0.3 dex precision (or perhaps better) for low-resolution spectra 
having S/N > 20. For all stars with multiple measurements of [Fe/H], we adopted the weighted mean from 
all measurements. Comparisons of [Fe/H] for the stars with multiple observations show a standard deviation 
of ~ 0.3 dex in the residuals between individual observations, suggesting that this is a reasonable value for 
the typical metallicity uncertainty. 



3. Virgo Substructure Candidate Selection 



We have already shown in Section 2. 1 that the stellar excess of blue, faint stars in SA 103 is not 
due to expected Galactic halo populations. Our estimate that ~ 30% of the stars within the CMD box 
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surrounding the turnoff region should be Virgo members suggests that roughly 15-20 of the 53 spectroscopic 
targets within that box should be part of the Virgo substructure, with a handful of additional RGB and AGB 
candidates at brighter magnitudes. Figure |2] shows the CMD for all stars with measured proper motions in 
SA 103, highlighting stars having spectra with S/N > 10. Stars are color coded by radial velocity to search 
for clumping in Vgsr- From these data, there appears to be a higher concentration of high velocity stars (the 
lighter colored points) in the boxed region of the CMD. This is consistent with prior studies claiming that 
the VSS is comprised of high velocity stars. 

Figure [3] shows histograms of the line-of-sight VgsrEI of all stars in the SA 103 field; the left panel 
depicts all stars that have spectroscopic data with S/N > 10 (the color-coded points highlighted in Figure |2). 
In the right panel we show only a subset of faint, blue stars corresponding to the boxed turnoff region in 
Figures [Hand |2] In the absence of substructure, stars within this box will be primarily Galactic halo turnoff 
stars. These halo stars should have a velocity distribution centered at Vgsr ~ km s -1 . A Gaussian fit to all 
halo stars from the Besancon model for the SA 103 field yields a velocity dispersion for the halo of ~ 90 
km s -1 , which we overplot in the right panel as a reference. The histograms in each panel are compared 
to the aforementioned Besancon model (the dashed blue lines) corresponding to the same field. In the left 
panel, the model was normalized so that the number of stars in the velocity range -300 < Vgsr < 100 km 
s" 1 is equal to the number of stars in the same range in SA 103. This velocity range was chosen so that the 
normalization includes only the region of the SA 103 data that is free from apparent substructure. In the 
right panel, the Besancon model counts were scaled by the fraction of the total number of stars in SA 103 
within the CMD-selection box that were observed; i.e., because we obtained spectra of 53 of 183 total stars 
in the turnoff box (determined from SDSS DR8 data), the Besancon histogram was scaled down by a factor 
of 53/183. If the SA 103 data contained only halo stars, the Besancon model scaled in this manner should 
fairly accurately predict the total number of stars in the right panel of Figure [3l clearly there are more stars 
in SA 103 than predicted by the model. There seems to be an excess of stars in the region of Vgsr ^ 100 km 
s -1 relative to both the Besancon model and the theoretical halo distribution. This exces s appears to have a 



narrow peak centered at a value somewhat greater than the (Vgsr) = 137 km s" 1 found by Brink et al. ( 20101) . 
though there also appear to be extra stars over a much bro ader range of high vel o cities. Note also th at the 
peaks at Vgsr ~ - 75 km s -1 and ~ -170 km s -1 seen by LNewberg et all (120071) . iBrink et all (120101) . and 



Casey et al. 



(12012 ) are apparent in this field of view; these peaks have been suggested to contain Sagittarius 



tidal debris. 

To verify the existence of the apparent high-velocity excess, a comparison sample of all stars between 
0.15 < g — r < 0.5 and 18.5 < g < 20.5 that were observed spectroscopically by SEGUE from 265° < I < 
283° and 50° < b < 68° (i.e., a > 300 deg 2 area surrounding SA 103) was extracted from SDSS DR8. A 
histogram of velocities for these stars is seen in Figure @] Because of the prominence of the Sagittarius 
stream at negative Vgsr in this region, we show only Vgsr > km s" 1 in the figure. For comparison, we 
overlay (as a dashed, blue line) the Besancon distribution from the smaller SA 103 field, scaled so that the 
number of stars between < Vgsr < 100 km s -1 is equal to that in the SDSS histogram. There may be 



5 Velocity in the Galactic Standard of Rest, denned as Vgsr = Vheiio + 9.0cosZ?cos/ + 232.0cosfcsin/ + 7.0sin£>. 




Fig. 3. — The left panel shows the measured Vgsr for all stars with spectroscopic data in SA 103 as the 
solid (red) histogram. The right panel contains only stars within the CMD box outlined in Figures Q] and 
[2] The blue dashed line in both panels represents the velocities predicted by the Besancon model in the 
field of SA 103 within the same color and magnitude ranges as the SA 103 stars. On the left, the Besancon 
distribution has been scaled such that the area under the histogram within the region between -300 < Vgsr < 
100 km s -1 equals that of our SA 103 sample. On the right, the Besancon model is scaled by the fraction of 
the total number of stars within the CMD-selection box among SDSS DR8 data in this field of view that were 
observed; i.e., since we obtained spectra of 53 of 183 total stars, the Besancon histogram was scaled down 
by a factor of 53/183 to match the number of SA 103 stars observed spectroscopically. In the left panel, 
we include only stars between 16 < g < 20,0.2 < g—r < 0.9 from the Besancon model to roughly mimic 
the color and magnitude range covered by the SA 103 Hydra spectra. On the right, a Gaussian centered 
at Vgsr = km s -1 with a a = 90 km s -1 (representing the typical Galactic halo distribution, with a = 90 
km s -1 derived from a fit to the stars identified as halo members in the Besancon model) is also superposed, 
again normalized to the area of the SA 103 histogram within the region mentioned above. In the right panel 
(i.e., within the main sequence turnoff region contained by the CMD selection box) there is a clear excess of 
SA 103 stars over both the model and the halo predictions at Vgsr 120 km s -1 . This excess is not confined 
to a narrow, kinematically-cold peak, but rather spreads over a large range in velocity. Note also that there 
may be stellar excesses present at -100 < Vgsr 5 ~50 km s" 1 and Vgs r ~ -170 km s" 1 ; these correspond to 
similar excesses seen in this region of the sky by lNewberg et all (120070 and lBrink et al.l (120101) and suggested 
to be related to the Sagittarius tidal stream. 



differences in the model between the smaller SA 103 field and the larger SDSS region, but such differences 
should be small because the expanded field is centered in SA 103. To represent expected Galactic halo 
populations, we also overplot a Gaussian centered at km s -1 with a a = 90kms _1 . As in our SA 103 data, 
the DR8 stars exhibit a clear excess in each histogram bin at high (Vgsr ^ 120 km s -1 ) velocities relative 
to the model and the predicted halo distribution. Because the SEGUE selection function has a patchy 
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Fig. 4. — Galactic rest-frame radial velocities of stars selected from 265° < / < 283° and 50° < b < 68° from 
SDSS DR8 within the same CMD box (i.e, between 0.15 < g-r < 0.5 and 18.5 < g < 20.5) as in SA 103. 
The Besancon model from the 1 deg 2 field of SA 103 is shown as a dashed blue line for comparison. Here, 
the Besancon model was normalized so that the areas of the histograms are equal within the region of 
< Vgsr < 100 km s -1 (the data was cut off at Vgsr > km s -1 because many Sagittarius stream stars 
appear in this region of sky at negative velocities). A Gaussian centered at km s -1 with a dispersion of 90 
km s -1 (representing the Milky Way halo) is plotted as a solid black line. There is a distinct excess of stars 
in every bin at Vgsr ^ 120 km s" 1 relative to both the Besancon predictions and the halo Gaussian. 

distribution in color and magnitude, these stars are likely not a complete representation of the velocity 
distribution in this region. The excess of stars at high Vgsr is present even in this sample, leading us to 
conclude that (a) the apparent excess in SA 103 is due to a real overdensity at high (Vgsr ^ 120 km s -1 ) 
velocities, and (b) the large spread of velocity in the excess is intrinsic, and not an artifact of the low 
resolution of our Hydra data. 

Once a rough idea of the velocity of Virgo candidate stars was known, we proceeded with a detailed 
candidate selection. The steps we undertook to identify a sample of Virgo substructure members were: 

1. Select all stars with 40 < Vgsr < 350 km s" 1 . This broad range in velocity was chosen to include all 
Virgo substructure members, but will also contain some MW contaminants. A total of 38 stars were 
selected in this manner. 

2. Narrow this sample to include only stars with at least 3 position measurements spanning at least a 40- 
yr time baseline to ensure that their proper motions are well-measured. In Figure [5] which compares 
proper motions of sample 1 in the left panel to those of the Besancon model on the right, there is 
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Fig. 5. — The left panel shows the proper motion vector point diagram (VPD) of all stars with well-measured 
PMs and 40 < Vgsr < 350 km s -1 in SA 103. On the right we show stars from five realizations of the 
Besancon model within the same range of Vgsr as the left panel. Color and magnitude cuts from 16 < go < 
20 and 0.25 < (g-r)o < 0.7 were applied to the Besancon points since spectroscopic data from SA 103 are 
roughly limited to within this region. The SA 103 stars predominantly occupy a much different region than 
the predicted Galactic halo stars in this field. The mean proper motion of the stars in the velocity-selected 
SA 103 sample is (/j, a cos5,fis) ~ (-5.2,0.0) mas yr -1 . 



a clear clumping of the proper motions about the median value of (fi a cos 5, /is) = (-4.8,-0.2) mas 
yr -1 . Furthermore, this clump is clearly offset from the expected PMs of halo stars predicted by the 
Besancon model. We thus choose to keep only those stars with Afi < 8 mas yr -1 , where Afi is the 
offset from the median PMs of the stars in selection 1 . This selection retains 30 stars. 



3. Calculate the offset of each star in sample 2 from an appropriate color-magnitude sequence for the 
samp le 2 stars . For this purpose we chose to use the empirical SDSS globular cluster ridgelines 
from 



An et al. 



(200 8J). This mi tigates any uncertainties that would arise from the use of theoretical 
isochrones (see lAn et all 120081 for illustration of the differences between empirical and theoretical 
ridgelines for MW clusters). Implicit in our use of globular cluster ridgelines is the assumption that 
the Virgo stars come from a predominantly old (> 10 Gyr) stellar population; this is a reasonable 
assumption if the stars came from a dSph galaxy. The median metallicity of the 23 stars from sample 
2 with spectra having S / N > 15 is [Fe/H] « -1.1 with <7 [Fe / H] = 0.6. We used the SDSS (g, g-r) 
ridgeline from lAn et all ( 2008 ) for the globular cluster with metallicity nearest this value - that of 
M5, a cluster with [Fe/H] = -1.27. We plotted the dereddened CMD of all stars from SDSS DR8 in 
the SA 103 field of view (the stars that make up the upper-right panel of Figure [D, highlighting the 
candidates in sample 2. We then overplotted the M5 ridgeline, shifting it in distance until we visually 
determined a "best-fit" distance of ~ 14 kpc. We excised all stars more than 0.2 magnitudes from 
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the 14-kpc M5 locus in (g,g — r) at magnitudes fainter than g = 17.5; this selection criterion serves 
to remove stars well off the expected CMD locus based on the observed turnoff, and the magnitude 
limit means that this selection is only applied for the turnoff, subgiant branch, and lower RGB. The 
upper RGB and AGB are not included because these regions are much more sensitive to differences in 
the (unknown quantities) age and [Fe/H]. Although our derived metallicity is higher than most other 
estimates for Virgo-substructure stars, and the distance we have found is slightly nearer than most 
findings, this has little effect on the CMD filtering. A generous width of ±0.2 magnitudes from the 
ridgeline would likely encompass all of the same stars if we were to change the parameters of the 
ridgeline used for filtering. This additional cut reduces the sample to 25 stars. 

4. After selection number 3, the proper motions of most of the remaining stars were clumped tightly 
about (ji a cos 5, us) = (-5.1,-0.2), with some obvious outliers at the edges of the distribution. This 
is illustrated in Figure [6j which shows the proper motion vector point diagrams (VPDs) of all stars 
with 40 < Vgsr < 350 km s -1 (i.e., sample 1) in the left panel, color-coded by their measured radial 
velocities. We noted that many of these stars in the outer regions of the PM distribution were also at 
the extreme ends of the Vgsr selection range, suggesting that they were Galactic interlopers, or had 
suspect proper motion measurements (typically having only 4-5 of the 14 possible position measure- 
ments). We thus chose to remove all stars with Afi < 5 mas yr -1 from sample 3, where the 5 mas 
yr -1 limit was chosen by examination of the left panel of Figure [6] The VPD of the remaining sample 
of 16 stars is shown in the right panel of Figure [6] None of the weighted mean proper motions and 
radial velocities of samples 3 and 4 varied more than la from the sample 2 measurements, suggesting 
that the somewhat arbitrary choices of exclusion criteria were not subjectively biasing our selections 
strongly. This final sample of 16 stars is highlighted in the CMD in Figure|7] The CMD shows all stars 
in sample 3 as filled symbols, color-coded by radial velocity. The 16 Virgo candidates constituting 
sample 4 are filled circles, and stars exclu ded by the proper motion cu t are triangles. The RR Lyrae 
star RR 167, a known VSS member from ICasetti-Dinescu et al.l (120090 . is the cyan-colored circle at 
g ~ 17, g-r ~ 0.2, and the t wo stars with Vgsr > 300 km s -1 are shown as green filled hexagons. 
The overlaid isochrones from lGirardi et al. ( 2004 ) are for 10 Gyr populations, with [Fe/H]=-1.2, and 
distances of 14, 17, and 20 kpc (spanning the range of most Virgo substructure detections from the 
literature). The 14 kpc isochrone appears to be the best match, and the 20 kpc ridgeline is clearly 
discrepant with the Virgo candidates. 



We note that the final sample of 16 stars contains three stars at g ~ 17, (g-r) < 0.4 that, if they are 
members, must be horizontal bra nch stars (specifically, RR Ly rae). One of these three stars, known as RR 
167, has already been shown by ICasetti-Dinescu et al.1 (|2009h to be a member of the VSS, so we believe 
that it is reasonable to include the other two in ou r sample. (Note th at although these stars may not have 
been detected as RR Lyrae in the QUEST catalog of IVivas et all (120041) . the ~ 70% completeness of QUEST 
for RRc-type stars in this region of sky suggests that some bona fide RR Lyrae may have been missed by 
QUEST.) There are also two stars in the sample at rather high RVs of Vgsr > 300 km s" 1 , shown as filled 
green hexagons at g = 19.75 and g = 19.91 ; these two stars have proper motions consistent with membership, 
and because there is not a narrow peak in the RVs due to Virgo stars, we choose to retain them in our samples. 
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Fig. 6. — The panels show proper motions of Virgo substructure candidates color coded by Vgsr- The panels 
on the left show stars with 40 < Vgsr < 350 km s" 1 , while the ones on the right show just the final selection. 
Again, RR 167 is highlighted in cyan. 



However, these two high-velocity stars appeal - to lie below the MSTO of the 14-kpc isochrone, and thus may 
be more distant than the majority of the Virgo members. This may signify extension of the Virgo structure 
along the line of sight, including possibly a velocity gradient with distance. 

Of the 16 stars in the "final" sample, 11 had spectra with S/N > 20, and all 16 had S/N > 14. From 
these, we derived metallicities with ~0.3-dex precision; Figure [8] shows metallicities for all SA 103 stars 
with spectroscopic S/N > 15 as a. function of velocity. The first selection of 40 < Vgsr < 350 km s -1 Virgo 
candidates are highlighted with larger (red/green) symbols; large red circles are the final Virgo sample, 
the two green circles are the stars flagged in Figure |2]for having Vgsr > 300 km s -1 , and the downward- 
facing triangles are stars removed from the final sample based on proper motion and CMD criteria. One 
star from the final sample does not appear in the figure because of what is likely a spurious [Fe/H]=-3.75. 
The 15 stars at [Fe/H]> -3.5 yield a mean value of ([Fe/H]) = —1.1 ±0.2, with (7[ Fe /jj^ = 0.7 (note that this 
becomes ([Fe/H]) = -1.0 ±0.1, o^ Fe j H ^ = 0.6 if the single star at [Fe/H]~-2.8 is excluded). This metallicity 
for the Virgo candidates is consistent with that found from isochrone fitting, and more metal-poor than the 
predominantly thick-disk stars toward the left side (i.e., with -150 < Vgsr < 50 km s -1 ) of Figure [8] 



4. Kinematics of Virgo Debris 



Statistical kinematical information for the 16 identified Virgo substructure candidates is given in Ta- 
ble |2] Shown in the table are the mean spatial coordinates and error-weighted mean proper motions, 
radial velocities, and metallicities. The best fit kinematical parameters were converted from proper mo- 
tion and radial velocity to r ight-handed Galactocentric Cartesian velocities using the method described in 



Johnson & Soderbloml (11987b . A Sun-Galactic center distance of 8.0 kpc was used, as well as solar velocity 
components in the Galactic rest frame of (U,V,W) = (9,232,7) km s -1 . The resulting Galactic Cartesian 
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Fig. 7. — Dereddened color magnitude diagram of all stars with 40 < Vgsr < 350 km s" 1 , -12 < \i a cos 5 < 2 
mas yr -1 , -6 < < 6 mas yr -1 , and measured on at least 3 photographic plates highlighted as circles and 
triangles. Symbols are color coded by the measured Vgsr- The greyscale Hess diagram in the background 
represents stars detected within the field usin g SPSS DR8 data. Is ochrones of age 10 Gyr, metallicity [Fe/H] 
= -1.2, and distances 14, 17, and 20 kpc from Girardi et al. ( 2004b are superposed in green, blue, and yellow, 
respectively. Stars that were removed from the sample based on CMD position and proper motion are shown 
as triangles. The s tars in the final selection ar e shown as filled circles. RR 167, an RR Lyrae star shown to be 
a VSS member (|Casetti-Dinescu et al.ll2009l) . is highlighted in cyan, and the two stars at high (Vgsr > 300 
km s -1 ) velocities are filled green symbols. 



velocity components of the Virgo debris are (U,V,W) = (-264.1,-331.1, 123.9) ±(65.7,50.5,33.0) km s _1 . 



These best fit kinematic pa rameters we r e use d to generate an orbital model. The Galactic gra vitational 
potential adopted is that from IWillett et al.l (120091) (which is, in turn, based on Law et all 12005): a three 
component potential that consists of a Miyamoto & Nagail (119751) disk, a Hernquist spheroid re presentin 
the bulge, and a logarithmic halo. The Galactic potential parameters are the same as those in Law 



ntmg 
et all 
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Fig. 8. — Metallicity ([Fe/H]) vs. Vgsr f° r ai l stars i n our sample having spectra with S/N > 15. The 
majority of the stars between -150 < Vgsr < 50 km s -1 are likely MW thin/thick-disk stars with [Fe/H]> 
—1.0. The large filled circles and triangles are the stars in our initial Virgo-candidate velocity selection. 
Triangles are stars that were removed from the sample by subsequent selection criteria, leaving only the 
large filled circles as the final Virgo candidates. The two points at the far right side have been colored 
green to distinguish these Vgsr > 300 km s -1 stars from the rest of the sample (as in Figure |2). The mean 
metallicity of the final sample is [Fe/H] = -1.1. 



(2005). Orbits were constructed using the mkorbit tool of the NEMO Stellar Dynamics Toolbox (ITeuben 
19951) . The orbit was integrated for 3 Gyr both forward and backward using the values tabulated. 
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Fig. 9. — The orbit resulting from our measured Virgo kinematics, integrated for 1 Gyr forward and back- 
ward, is shown as the thick blue line in each of the panels above. The dashed portion is the backward 
integration (i.e., the path of the infalling progenitor), and the solid part represents the forward integration of 
the orbit. The left panel depicts the orbital path in equatorial coordinates, and the right panel in Galactic co- 
ordinates. The kinematics used to generate this orbit are given in Table [2] The location of SA 103 is shown 
as a large black dot at (a, S) = (178.7,-0.6)°, or (l,b) = (274.6,59.1)°. The positions of known globular 
clusters and dwarf galaxies are plotted in the background as yellow and green points respectively. 

4.1. Orbit 

Figure [9] shows the resulting orbit in both equatorial and Galactic coordinates. Milky Way dwarf 
spheroidal galaxies (dSphs) and globular clusters (GCs)0 are overplotted for later reference with respect 
to finding possible orbital associations and/or the progenitor. The orbit has apo- and peri-centers of R a = 
52^4 kpc and/? p = 5.6+Qg kpc, with eccentricity e = (R a -R p ) /(R a +R p ) = 0.81^ 27 (errors on these quantities 
are derived by propagating the uncertainties in the kinematics). The orbital inclination with respect to the 
Galactic plane is i = arcsin(|z max |//? zmax ) = 54 ±3°, where |z max | is the maximum height above or below the 
plane, and R zmax is the distance from the Galactic center where the orbit passes through |z max |- The orbital 
period (measured between successive apocenters) is P OI \, = 0.7+Q3 Gyr. Plots of the orbit in distance and 
Vgsr vs. Galactic longitude are shown in Figure [10] Again the dwarf galaxies and globular clusters are 
plotted in the background for reference. The high eccentricity (e ~ 0.8) of this orbit is evident, with the 
satellite's trajectory having just passed pericenter after swooping in from the outer halo. 

The orbit was converted to Galactic right-handed Cartesian XYZ coordinates, with the convention that 
the Sun is centered at {X,Y,Z) = (-8,0,0) kpc and the direction of the X-axis points from the Sun towards the 
Galactic center. Plots of the derived orbit projected onto the Galactic XY and XZ planes are shown as solid 
(blue) lines in Figure [TT] In these projections, it is clear that the observed Virgo structure has just passed 
perigalacticon in its orbit after swooping in toward the Galactic center from a distance of > 50 kpc on a rather 
radial orbit. To give a sense of possible orbital paths of debris due to both the uncertainties in our orbital 



s From the compilation of Dellomo & Newbera l 20091) . available at jhttp : //homepages ■ rpi . edu/~newbeh/mwstructure/MilkyWay£ 
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Fig. 10. — Derived orbit as a function of distance and Vgsr- The upper panel shows heliocentric distance 
vs. Galactic longitude for the derived Virgo substructure orbit (thick blue line). In the lower panel, we give 
Vgsr vs. Galactic longitude. Plots are labeled similarly as in Figure |9j with Galactic satellites (dSphs and 
GCs) shown as green/yellow dots. It is evident that the current Virgo debris is near orbital perigalacticon 
and has just passed fairly near the Galactic center. 

parameters and the dispersion of 3-D velocities in the VOD, we show 10 orbits generated by randomly 
picking from within la Gaussian distributions of the uncertainties on each of the orbital parameters in 
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Table[2j Most of these "perturbed" orbits do not depart drastically from the mean orbit. 



These plots are similar in appearance to the orbits shown in ICD09l Our orbit, derived from 16 stars 
rather than a single RR Lyrae, has similar eccentricity and apogalacticon as i n the pr evious work, confirming 
that we have measured kinematics of debris from the same structure as the ICD09I study. Note also that on 
first glance our data suppor t (within the uncertainties) an association of NGC 2419 with the Virgo structure, 
as was suggested by |CD09| . This and other plausible orbital associations will be explored in the following 
section. 



4.2. Possible Orbital Associations 



To search quantitatively for possible orbital associations of satellites to our derived orbit, we compared 
the locations and velocities of known GCs and dSphs with the path and velocity of the mean orbit. This was 
done by defining a chi-squared value made up of the sum of the residuals (scaled by the uncertainties) in 
position and velocity between each point of the orbit and each Milky Way satellite's known Galactic XYZ 
position and Vgsr velocity. The uncertainties in spatial positions at each orbital time-step were taken to be 
30%, and in Vgsr to be 30 km s -1 . These large ranges should encompass the uncertainties in the kinematics 
used to generate the mean orbit0 The minimum x 2 is located for each object, with the assumption being 
that this represents the nearest point in phase space to the satellite. Globular clusters within 15 kpc of the 
Galactic center were excluded from the analysis because the Virgo substructure's orbit will pass relatively 
near all of the globular clusters that are located within this region. Without 3-D kinematics for all of these 
inner GCs, it is difficult to assess their possible kinematical association with the Virgo orbit. 



The globular cluster NGC 2419, which was claimed by |CD09| to be plausibly associated with the Virgo 
orbit, is located ~ 30 kpc from the apocenter of our orbit. Note, however, that a slightly larger space 
velocity for Virgo debris at the SA 103 position would send an orbit through almost exactly the position of 
NGC 2419 (see, e.g., the maroon-colored "perturbed" orbit in Figure [TTb. The association of Virgo debris 
with this globular cluster is thus tentatively confirmed by our data. 

The lowest \ 2 from our search (and, in fact, the o nly match with y 2 < 5.0) was for "SEGU E lb", which 
denotes the apparent velocity substructure found by iGeha et al.1 (120091) and lSimon et all (1201 lb in the field of 
view of the ultra-faint dSph SEGUE 1. SEGUE 1 itself is spatially a good match to our Virgo substructure 
orbit, but the radial velocity predicted by our orbit agree s much better with the debris at ~ 100 km s -1 higher 
velocities (and roughly the same distance, according to iSimon et al.ll2011[) than with SEGUE 1 itself. 



In this section, we briefly discuss possible satellite associations found by our fitting technique as well as 
other known stellar overdensities; we will revisit each of these particular cases one-by-one in the following 



The assumption of a constant fractional uncertainty in the distance and fixed velocity uncertainty at all positions is chosen to 
simplify these calculations. In fact, for a given set of initial conditions and their associated uncertainties, the position and velocity 
uncertainties vary substantially at different points along the path of the orbit. Simulating this accurately would require a more 
sophisticated approach that is beyond the goals of this simple search for possible orbital associations. 
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Fig. 1 1 . — The upper panel shows the projection of the Virgo structure orbit onto the Galactic XY plane, and 
the lower panel gives the XZ projection. Thick blue lines are the orbit derived from our mean kinematics; 
solid blue lines are the forward integration, and dashed blue lines the backward integration. The thinner lines 
of other colors represent orbits generated by selecting from a Gaussian I — a distribution of the errors in the 
three velocity components and the distance, and illustrate the constraints placed on the Virgo progenitor 
orbit by our measurements. The Sun is at (X,Y,Z) = (-8,0,0) kpc in this coordinate system, with positive X 
pointing in the direction from the Sun to the Galactic center. The position of SA 103 is represented as a large 
black circle. We also highli ght as a large black dot the globular cluster NG C 2419, which was te ntatively 
associated with the VSS by ICD09l along with the "90 kpc Sgr debris" from INewberg et all (120031) as plus 
signs. Two other substructures found to lie near our derived orbit - SEGUE 1 (or, more specifically, debris 
in the field of view of this ultra-faint dSph) and the Pisces Overdensity - are also given as large black dots. 
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section after we present the Af-body model based on our orbit. The two main structures that were found to 
be consistent with the integrated orbit were NGC 2419 and the velocity substructure in the field of SEGUE 1 
(which we denote "SEGUE lb"). Another structure that may be related to the Virgo structure that is neither 
a GC nor a dSph is the Pisces Overdensity. This stellar overdensity is apparently located near a turnaround 
point in its orbit, and possibly contains distinct moving groups at positive and negative radial velocity, which 
further suggests that it is debris near apocenter of its orbit. An additional possibly-related structure is the 

to be part of Sagittarius located at 90 kpc from the Galactic 
center, in close proximity to NGC 2419. These structures are briefly introduced here, and will be discussed 
in more detail in the following section of this work. 



debris previously claimed (Newbers et al 



4.2.1. NGC 2419/ 90 kpc Sagittarius debris 



NGC 2419 is located at (l,b) = (180°, 25°) with heliocentric and Galactocentric distances of 84 and 
91 kpc respectively. Its helioc entric radial velocity i s -20 km s" 1 , corres ponding to a Galactic standard 
of rest velocity of -14 km s -1 (|Peterson et al.lll986T) . iNewberg et al. d2003h find an excess of A-type stars 
extending at least 10 degrees around NGC 2419 at about the same quoted distance. The authors claim 
that since NGC 2419 is contained within the debris, the tidal debris can be inferred to be at apogalacticon 
of its orbit, consistent with their findings of being coincident with a portion of the Sgr stream near its 
apogalacticon. However, this was merely a tentat ive association of these stars with Sagittarius; our current 
understanding of the Sgr tidal debris system (e.g.. Law & Majewskill2010i) suggests that there should not be 
debris at 90 kpc. 



4.2.2. SEGUE 1 



First discovered by iBelokurov et all (12007b using SDSS imaging data, this dwarf galaxy has subse- 
quently been the subject of close examination because of its peculiar properties, most notably its appar- 
ent high mass-to-light ratio. Located at (RA,Dec)= (152°, 16°) and (l,b) = (220.5°, 50.4°) at a distance 
of cfhei = 23 ±2 kpc and dec = 28 kpc, it is a good match to our derived orbit of the Virgo substructure. 
Belokurov et al. describe an apparent tidal tail feature extending from the northeast to the southwest in 
eq uatorial coo rdinates. However, subsequent studies have not been able to reproduce such results. A study 



by lGeha et al. 



EIMOS spectroscopy presented radial velocity and metallicity data, and 
found a main grouping of stars at about Vqsr = 111 ± 1.3 km s -1 with a sma ller group of 4 stars at ~ 200 
km s . They also obtain [Fe/H] = -3.3 ± 0.2 for one of the SEGUE 1 stars. ISimon et al.l (|201 1|) extended 
upon this survey using mor e data and obtain similar numbers to those of Geha et al. They also confirm the 
~ 200 km s" 1 feature from beha et all (120091) . identifying 24 stars within the peak at (Vgsr) = 203.8 ± 1.7 



km s l , compared to the original 4 s t ars. T hese 24 stars represent a kinematically cold population, with 



cry = 7.0 ± 1.4 km s l . ISimon et all (1201 11) suggest that the stars making up this feature must be more 



metal-rich than those of SEGUE 1 . 
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Our initial search for positional and kinematical associations to our Virgo orbit found SEGUE 1 was a 
fairly good match. However, we noticed that the ~ 200 km s" 1 feature actually fits in better with our findings 
at the given location. Thus the chi-square test was rerun with an additional element, labeled "SEGUE lb", 
to distinguish it from the main SEGUE 1 population. The SEGUE lb stars are an excellent match to our 
orbit in all three spatial dimensions (within 1.3 kpc), and offset from the best matching point by only ~ 10 
km s -1 in radial velocity. 



4.2.3. The Pisces Overdensity 



The Pisces Overdensity was first found by ISesar et al.l (12007b as an overdens ity (dubbed "S t ructur e J") 



Watkins et al. 



(2009) and 



of distant RR Lyrae stars in SDSS Stripe 82. It was subsequently examined by 
shown to stretch across an area from 63° < / < 93° and -60° < b < -46° at a Gala ctocentric distance of 
79.4 ± 14.1 kpc and a metallicity of [Fe/H]= -1.48 ±0.28. iKollmeier et all (120090 measured the radial 
velocities of RR Lyraes within SDSS Stripe 82 and found two distinct groups of stars, one centered at a 
heliocentric radial velocity of -75 km s" 1 and the oth er around -198 < V r < -155 km s" 1 . They find similar 
distances and metallicities as Watkins et al. Later, ISesar et al. I (hoiObh expanded the selection of stars of 
Kollmeier et al. and confirmed two distinct moving groups of Galactocentric velocity of -52 ± 11 km s -1 
and 50 ±3 km s . This finding of two separate velocity groups in opposite directions suggests that this 
overdensi ty is at a turn i ng po int in its orbit. This too agrees nicely with the orbit that we have obtained in 
this area. ISharma et al.l (120101) find a structure among 2MASS-selected M giants in the vicinity of the quoted 
location of the Pisces Overdensity, which is labeled "A16" in their paper. They find that the peak density lies 
just outside the range of SDSS Stripe 82. They hypothesize that the Pisces structure is an extended remnant 
of a satellite that was on a highly eccentric orbit and got tidally disrupted. However, in the same paper, they 
were not able to detect VSS members, calling into question such an association with VSS. 



5. Discussion: the Progenitor of the Virgo Stellar Substructure 

5.1. N-body Simulation 

With the orbital parameters of Virgo debris well-constrained, we turn to the question of what type of 
progenitor could produce a structure with properties like those of the Virgo substructure. With the extant 
spatial and kinematical data, it is unclear whether all of the numerous detections of substructure in Virgo are 
related. However, we will show that a simple model based on our precise kinematics explains nearly all of 
the observed properties. 

To understand the orbital dynamics of the Virgo substructure progenitor more clearly, semi-analytical 
Af-body simulations were generated using our derived orbit. The kinematical parameters determined for 
Virgo debris in SA 103 were used to generate orbital models using the same gravitational potential as the 
test particle orbit. The simulation templates used were 10,000 particle Plummer spheres placed on our 
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derived orbit for the Virgo structure. The satellites were evo lved in the Ga lactic potential for 3-4 Gyr using 
the gyrfalcON tool of the NEMO Stellar Dynamics Toolbox ( Dehnen|[2002 ). and their ending time tailored to 



place the core of the progenitor near the current VOD. This strategy of placing the progenitor (or its remains) 
at the position of the VOD is motivated by our derived orbit. Pile-ups of tidal debris are typically located at 
the apo- or peri-center of a satellite's orbit; however, we showed in the previous section that our measured 
kinematics place Virgo substructure stars at an orbital phase just past pericenter. Thus, in order to have a 
cloud-like pile-up of stars at the VOD position (i.e., at an orbital phase that is not at peri- or apo-center), the 
stellar overdensity is likely the remains of the progenitor itself. 

Because the progenitor of the Virgo structure is unknown (if it even still exists), we must explore the 
behavior of disrupting satellites of different masses. Initially, we placed a globular-cluster-sized object with 
M = 1 x 10 5 M Q and Plummer radius r p i = 50 pc on our derived orbit. On the eccentric orbit derived for 
Virgo stars, the globular cluster is rapidly shredded upon its pericentric passage. The A^-body results of 
this globular cluster progenitor are seen in Figure [12] which focuses on the region of sky around SA 103. 
The satellite is already extended into a narrow tidal stream after its first perigalactic passage (in the upper 
panel), and is completely disrupted by around 3 to 4 Gyr, corresponding to the third perigalactic passage. 
This creates a narrow, well defined stream in the sky, contrary to observations of a large cloud-like feature 
seen in Virgo. We note also that remnants of disrupted globular clusters are typically rather kinematically 
cold; the absence of a narrow velocity peak among the Virgo substructure data thus also argues against a 
globular cluster progenitor. Therefore, a globular cluster-like object is unlikely to be the progenitor of the 
substructure(s) in Virgo. 

Having ruled out a globular cluster as the Virgo substructure progenitor, we then tried a larger, more 
massive dwarf galaxy -like satellite. Specifically, we first simulated a dSph of mass similar to that derived by 



Newberg et all (120101) for the Orphan stream progenitor (M = 3 x 10 6 M Q ), configured in a Plummer sphere 
of radius r v \ = 0.2 kpc. This low-mass dwarf galaxy was run through an A-body simulation with the same 
MW parameters as the globular cluster model, on our derived Virgo substructure orbit. As with the globular 
cluster, the Orphan-like progenitor disrupts quickly on this destructive orbit. As can be seen in Figure [13] 
this comparatively kinematically "warmer" object dissolves much more rapidly than the globular cluster. 
The remnant becomes more "puffy" than the disrupted globular cluster, but still remains embedded in a 
narrow stream. This more closely resembles what is seen in Virgo, and suggests that an even larger, more 
massive dSph may produce a large-area, cloud-like overdensity when disrupted on this orbit. 

Finally, we simulated a 10,000 particle, Sagittarius-sized template dwarf galaxy (a Plummer sphere 
with mass of M ~ 10 9 M Q and scale radius r p i = 0.9 kpc) on the Virgo test-particle orbit. The results 
of this gravitational A-body simulation are displayed in Figures [14] (spatial distribution), [15] (Vqsr and 
distance vs. longitude), and [16] (polar projection of the spatial positions). The disruption of a Sgr-sized 
object on this orbit produces several structures that span large (hundreds of square degree) areas of the 
sky (Figure [14). The large, diffuse overdensities centered at (l,b) ~ (60°, -50°) [(a, S) ~ (330°, -5°)] and 
(J,b) ~ (210°,40°)[(a,<5) ~ (110°, 20°)] are the results of tidal debris piling up at apogalacticon. These 
distant debris (at Rqc ~ 80 and ~ 50 kpc, respectively; see Figure \\5} piling up at the turning points of the 
disrupting dSph's orbit are expected for an object on such a radial, destructive orbit (see, e.g., the "umbrella- 
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Fig. 12. — Shown here is an Af-body simulation of a globular cluster sized object evolved along the derived 
test particle orbit, in 4 different time slices, spaced 1 Gyr apart. The cluster is modeled as a Plummer sphere 
with mass of M = 1 x 10 5 M Q and Plummer scale radius of 50 pc, represented by 10,000 red particles of 
equal mass. The object is seen to rapidly stretch out into a thin filament, falling apart between 3-4 Gyr, or 
after 3 perigalactic passages. 



like" structures in lJohnston et al.H2008h . One may expect a dSph seen shortly after a perigalactic passage to 
have a distorted morphology with well-defined tidal tails emanating from it. However, the A-body model 
guided by our precise kinematics of Virgo stars shows that the large tidal shock imparted to a Sgr-like 
progenitor upon its passage < 10 kpc from the Galactic center causes the dSph to inflate to an apparent size 
(at d = 15 kpc) of hundreds of square degrees (note that we are not concerning ourselves about whether the 
particles are bound to the progenitor at this stage; most of the stars must be unbound at this point in the 
object's evolution). This is seen more clearly in the left panel of Figure [16] which shows a polar plot of 
Af-body model debris in Galactic coordinates, centered on the nort h Galactic po l e. Giv en for comparison in 



the right panel of Figure [16] is a reproduction of Figure 5 from iNewberg et all (120071) . in which the Virgo 
substructure is clearly seen as an overdensity of photometrically-selected F-turnoff stars from SDSS. The 
progenitor (or its remains) in our model spans roughly the same area as the overdensity seen by Newberg et 
al., centered on the same position, including possibly an elongated shape along Galactic latitude (i.e., from 
left to right in the figure). 

This structure, seen at roughly the position of SA 103 (denoted by the large black dot in Figure [14] 
at (l,b) ~ (275°, 59°), or (a, 5) ~ (179°, -1°)) and spanning ~ 500 square degrees, contains ~ 16% of 
the original 10,000 particles making up the progenitor. Since the progenitor had a mass of 10 9 M Q , this 
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Fig. 13. — Results of an Af-body simulation of an Orphan-like progenitor (a dwarf-galaxy-sized object) 
evolved along our derived Virgo structure orbit. The panels show snapshots at the position of SA 103 
in four time slices spaced ~1 Gyr apart. This object was configured as a Plummer sphere with mass of 
M = 3 x 10 6 M Q and a scale radius of 200 pc. The satellite dissolves quickly (in less than 3 Gyr), with the 
remnant showing a narrow stream-like structure at early times and virtually no spatial concentration at later 
stages. 



corresponds to a mass of ~ 1.6 x 10 8 M© in the area surrounding the remains of the progenitor. As an 
exercise, we estimate the luminosity corresp onding to this mass by assuming a mass-to-light ratio typi cal of 



Milky Way dSphs@ (f ) = 100 (see, e.g.. Iwolf et aDboiol . Iwalker et alJl2009l Isimon & Gehall2007h : this 



yields a V-band absolute magnitude of M v = -10.7 (or M v = -13.2,-8.2, respective ly, for f = 10, 1 000). 



This compares favorably to the estimate of the total integrated light of the VOD bv Prior et a l. (2009 ) in a 



760 d eg z area of My = -1 1 .9 (or My = -10. 1 using the more nearby lJuric et al.H2008l distances). lJuric et al. 
(|2008l) found a total r-band luminosity of 0. 1 x 10 6 L Q , or M r = -7.8, for the overdensity assuming a distance 
of 10 kpc and 1000 deg 2 area. Thus, the remains of the Virgo progenitor as determined by our derived orbit 
are consistent with the observed stellar overdensity that covers hundreds of square degrees. 

Another mystery surrounding the numerous overdensities detected in Virgo is the disagreement among 
radial velocities between studies of debris at different positions. As we showed in Section|3l even in SA 103, 
which spans only 40' x 40', there is a significant excess of high-VbsR stars that forms a broad peak (see 



In our model, the object is probably no longer a bound system, making the assumptions of these calculations dubious. We 
present these results simply as an attempt to compare A'-body particles to the observed stellar overdensity. 
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Fig. 14. — The results of an A^-body simulation of a Sagittarius-sized (M = 1 x lO 9 M0,r p i = 0.9 kpc) Plum- 
mer sphere on our derived Virgo substructure orbit. The left panel depicts the debris in equatorial coordi- 
nates and the right panel shows the simulation results in Galactic coordinates. The black point represents 
the present location of the observed Virgo debris in SA 103. This A-body simulation was initially placed 
3 Gyr in the past and integrated forwards for 2.95 Gyr. Three overdensities are immediately apparent: The 
ones at / = 60° and / = 210° are due to piling up of debris at apogalacticon. The third overdensity located at 
/ = 270° is progenitor related; in other words, the location of this cloud of debris changes depending on the 
duration of the simulation. We tailored the run-time of this particular simulation to place the progenitor (or 
its remains) at this position. 

Figure |3). This is more obvious among SDSS RVs shown in Figure 0]- a broad excess relative to the 
Besancon model is apparent at 100 < Vgsr < 250 km s -1 . This is much too broad to be the signature of 
a kinematically-cold tidal stream. If the VOD were due to the piling up of tidal debris at the apocenter of 
a satellite orbit, then a large velocity dispersion could be expected; however, the mean velocity at such a 
turnaround point must be Vgsr ~ km s -1 rather than ~ 150 km s" 1 as observed. Furthermore, we have 
shown that stars in the Virgo substructure are not at apogalacticon, having rather just passed their peri- 
galacticon. Examination of the upper panel of Figure [15] however, shows that our A-body model offers 
an explanation for the large spread in velocity of Virgo debris. The progenitor, having just experienced a 
close encounter with the Galactic center, has been kinematically heated so that its velocity dispersion has 
increased to ~ 20 km s -1 or more and its RVs span ~ 100 km s -1 . Note that the model seen in Figure [T5ldoes 
not predict debris with Vgsr > 200 km s -1 near the progenitor, as is seen in both the SA 103 data (Figure [3]) 
and SDSS (Figure HJ) as a long tail to high velocities. It is possible that additional substructure is present in 
this region of the sky. Alternatively, this long tail may be a feature of the disrupting Virgo progenitor that our 
simple model has not addressed - further modeling exploring the progenitor's mass and mass distribution 
may be able to reproduce this feature. 

To summarize, our Af-body model of a Sagittarius-sized dwarf galaxy with mean kinematics of the 16 
stars we determined to be Virgo substructure members qualitatively reproduces many of the observed prop- 
erties of the stellar excess in Virgo. We have shown that a large, cloud-like stellar structure like that in Virgo 
is likely to have originated in the tidal disruption of a massive dwarf galaxy. No dwarf galaxies are found 
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Fig. 15. — Af-body results from the simulation with a Sagittarius-like progenitor, as in Figure [14] The upper 
panel shows Vgsr vs Galactic longitude while the lower panel shows heliocentric distance vs. Galactic 
longitude. The model runtime was tailored to place the progenitor (and/or its remains) at the position of 
SA 103 (the large black dot). 

along the integrated path of the orbit we have derived, leading us to conclude that the Virgo substructure 
is itself the remains of the progenitor. The slight mismatch in the position of the disrupted progenitor in 
our model and that of the observed overdensity in Virgo is a consequence of the model runtime and our 
lack of knowledge of the central position of the stellar overdensity, which may lie outside the observed 
SDSS footprint. Details such as the size of the overdensity and the direction of its elongation may also be 
slightly discrepant between the model and observations, but are in rough qualitative agreement. Many of 
these differences are likely due to the fact that the initial conditions used to generate the model - namely, 
the kinematics of debris in SA 103 - may not actually correspond to the center of the progenitor. A more 
comprehensive exploration of parameter space (including, but not limited to, the mass of the progenitor, its 
position, and its initial density configuration) could likely resolve some of these issues, but is beyond the 
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Fig. 16. — Left panel: the Sagittarius-sized A^-body simulation in polar Galactic coordinates, centered at the 
North Galactic Cap. Here an overdense region around (l,b) = (285°, 55°) is clearly visible. Approximate 
boundaries of the possible associations are also drawn in the plot. Right panel: Density map of F-turnoff 
stars in the Northern Galactic Cap from Newberg et al. 2007. These stars were obtained from SDSS DR5 
based on the criteria 0.2 < (g— r)o < 0.3, (w-g)o > 0.4, and 20.0 < go < 21.0. Darker areas contain more 
F turnoff stars. Comparing this figure to the left panel, one can see an apparent agreement between the 
observed location of the Virgo debris and that predicted by the N-body model. The Virgo overdensity is 
visible at the lower right edge of the data in the right panel. 

scope of this work. It is clear, however, that the morphology, spatial extent, luminosity, and large velocity 
dispersion of the observations are all qualitatively consistent with our simple model. We thus suggest that a 
massive dwarf galaxy was the progenitor of the substructure in Virgo. 

5.2. Other Possible Progenitors or Associated Overdensities 

In Section 14.11 we discussed a search for Milky Way satellites and stellar overdensities that may be 
associated with the Virgo substructure. While initially the search was conducted in hopes of identifying 
the progenitor of the VOD, our modeling results suggest that the VOD itself is likely the remains of the 
progenitor. However, the infalling dwarf spheroidal that produced the VOD may have had satellites of its 
own, similar to the globular clusters associated with the Sgr and Fornax dSphs. If the progenitor had globular 
clusters orbiting in its potential, these should be found along the orbit of the progenitor. Also, since the VOD 
is near the pericenter of the progenitor's rather destructive orbit, there should be additional tidal debris that 
has collected at the apogalactica of the orbit, as seen in Figures Q4] and [15] In this subsection, we revisit 
the satellites and stellar substructures we identified in Section I4.2L and discuss them in the context of the 
A^-body model results. 
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Fig. 17. — Af-body model debris in the southern Galactic hemisphere. An overdensity of model debris is 
seen that is nearly coincident with the Pisces Overdensity between 60 < / < 90°. These stars, located at 
distances of ~ 80-90 kpc, correspond to the apogalacticon of our derived orbit. 

5.2.1. NGC 2419 and Nearby 90-kpc Stellar Overdensity 



The distant Milky Way globular cluster NGC 2419 was suggested by CD09 to be plausibly associated 
with their orbit of the VSS based on a single RR Lyrae star. The orbit passes near NGC 2419 at the next 
apocenter ahead of the Virgo substructure along its orbit. The cluster has a radial velocity of (Vgsr) = _ 14 
km s -1 , which, because it is close to zero, is consistent with NGC 2419 being at the turnaround point of its 
orbit. We find that our orbit does not reach the ~ 80 kpc distance of NGC 2419, instead turning around at 
~ 50 kpc at its nearest apocenter. However, within the uncertainties on the orbital parameters (as illustrated 
by the different orbits in Figure QT] with kinematics selected from within the 1 — a uncertainties), the orbit 
could plausibly be associated with a satellite or debris at such distances. NGC 2419 is an unusually large, 
luminous globular clust er, which has led to speculation t hat it is actually the stripped core of a former dwarf 
spher oidal galaxy (e.g., van den Bergh & Mackey 2004 ; Mackey & van den Bergh 2005 : Cohen et aDboiO . 
201 Vs . One might thus suggest NGC 2419 to be the progenitor of the Virgo substructure. We believe 
this is unlikely, because NGC 2419 has a metallicity of [Fe/H] = -2.1 (|suntzeff et al.lll988h . making it 
more metal-poor than the debris we have studied in SA 103, and even more deficient than the RR Lyrae 
studies have found for VOD/VSS stars. Because metallicity gradients in dwarf galaxies lead to more metal- 
rich populations being concentrated at their cores, tidal streams from disrupted d Sphs should have mean 
metallicities that decrease with distance along the stream from the progenitor (e.g. IChou et al.ll2007l) . It is 
hard to imagine how, if NGC 2419 is the progenitor, it could have a mean metallicity lower than that of 
all known d ebris. It is also possib le that the cluster fell in with the (dwarf-galaxy) progenitor. However, 



we note that 



PC 

(2( 



van den B ergh (2007) pointed out that the globular clusters associated with dwarf galaxies are 
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faint relative to the general globular cluster population; this would appear to rule out the extremely luminous 
NGC 2419 as a globular-cluster companion of an infalling dwarf galaxy. 



Furthermore, INewberg et all (120031) identified an overdensity of A-type (primarily BHB) stars at a 
mean distance of 83 kpc spanning over 200 deg 2 in the region surrounding NGC 2419. Very few BHB stars 
should be present that far out in the halo, and it was originally suggested by Newberg et al. that these stars 



are Sa gittarius tidal debris. In our A-body model of the Virgo progenitor, the position of the INewberg et al. 



(120031) stars on the sky is consistent with the collection of Virgo-structure stars at the apocenter of its orbit. 
We note that the A-body model produces debris about 40% closer than either NGC 2419 or the Newberg et 
al. stars, but this is not troubling. The orbit traverses distances consistent with both NGC 2419 and the A- 
star overdensity, and we have seen that A^-body models using slightly different kinematics (from within our 
error distribution) would produce model debris at those distances as well. We thus suggest that NGC 2419 
and the "cloud" of A stars nearby it may be associated with the Virgo progenitor. 



5.2.2. Substructure in SEGUE 1 Field of View 



Another of the close matches to our orbit is the dwarf galaxy SEGUE 1, which lies directly along 
the path of the Virgo structure orbit. However, we are not suggesting here that the extremely metal-poor, 
dark matter-dominated SEGUE 1 dSph is a s sociat ed with the Virgo structure. Instead, it is the "extra" 
radial velocity peak identified by iGeha et al.1 (|2009L with 4 stars), later s een more strongly by INorris et al. 
(|2010L 8 stars), and finally detected with 24 stars by ISimon et all (|201 lh that we believe to be tidal debris 
from the Virgo system. Because the spectroscopic targets were selected by Simon et al. to be consistent 
with the SEGUE 1 distance (23 ±2 k pc), the stars in the additional RV peak are presumably at roughly 
this same distance (|Simon et al.ll201ll quote a distance of ?a 22 kpc). However, their mean velocity is at 
Vgsr = 204 km s , offset by ~ 100 km s -1 from the SEGUE 1 velocity. Our A-body model predicts Virgo 
debris should be present at this position (l,b = 220,50°) at a distance of ~ 25 kpc (see Figures IT41 and [T3Ti. 
Furthermore, the upper panel of Figure [15] shows that the debris in our mode l is at 150 < Vqsr < 200 km 



s l , or almost exactly where the 24 "excess" stars in the ISimon et al.l (120111) study appear. We therefore 



conclude that the debris in the SEGUE 1 field of view i s plausibly a s sociat ed with the Virgo stream system. 
We note also that the mean metallicity is suggested by ISimon et all (|201 II) to be [Fe/H] = -1.3 for these 24 
stars in the "300 km s -1 stream". This is consistent with the metallicity we have measured for debris in SA 
103. Lending furt her credence to ou r claim of an association with the stream from the Virgo progenitor is 



Simon et all (|201 lh and Nie derste-Qstho l t et al . (2009) see an east-west extension of the 



the fact that both 

300 km s -1 stream stars. In fact, Figure 8b of ISimon et al.l (1201 11) shows the stars in this velocity feature to 



possibly be oriented along a southeast-to-northwest direction, which is exactly the orientation of the debris 
from our model of the Virgo progenitor. Based on all these facts taken together, we believe that the 300 km 
s" 1 stream in the SEGUE 1 field of view is debris associated with the Virgo substructure. 



-35- 



5.2.3. Pisces Overdensity 



Finally, we note that the "Pisces O verdensity", a s t ellar excess spanni ng a large area on the sky at a 
Galactocentric distance of 79 ± 14 kpc (ISesar et al.ll2007l : IWatkins et al.ll2009i) . is located at the position and 
distance of the most recent apocenter of our Virgo orbit. This overdensity was first found among RR Lyrae 
in SDSS Stripe 82 data at (a, 5) ~ (355°, 0°), or (l,b) ~ (88°, -58°), and subsequently seen to stretch to 
(l,b) ~ (105°, -53°) in 2MASS-selected M-giant stars (|Sharma et alJboioh . Examination of Figures [Hand 
[15] shows a pile-up of model debris at the position and near the distance roughly corresponding to the Pisces 
Overdensity. This is also illustrated in the polar view of model debris in the southern Galactic cap seen in 
Figure [T7] which overplots the approximate location of the Pisces Overdensity atop our model. Because 
this substructure has been detected in metal-poor RR L yrae (ISesar et al.ll2007l : IWatkins et aL||2009j as well 
as relatively metal-rich 2MASS-selected M-giant stars ( Sharma et al. 201ob . it has already been suggested 
to be related to a dwarf-galaxy disruption event. It is therefore remarkable that our /V-body model for a 
disrupting dSph on the Virgo orbit predicts just such a "cloud" of tidal debris as has been observed. The 
radial velocity one would expect at the turnaround point of the orbit is Vqsr = km s" 1 , with some large 
spread about zero velocity due to the viewing angle roughly along the path of both the outward-streaming 
and inward-bound debris (see Figure [l5]for an illustration of the large range of ex pected Vqsr in a ~ 20° 



range of Galactic longitude). Spectrosco py of RR Lyrae in th e Pisces structure by iKollmeier et all (|2009|) 
was combined with additional spectra by ISesar et al.l (|2010a|) to arrive at a sample of 12 RR Lyrae in the 
overdensity. These 12 stars separate into two disti nct peaks at (Vgsr) = 50±3 km s -1 with a = 10±3 km s -1 , 
and (Vqsr) = _ 52 ±11 km s -1 , a = 23 ± 9 km s -1 . Isesar et al ( 2010a ) considered the possibility that the two 
groups are drawn from the same Gaussian distribution, and concluded that the hypothesis of normality can 
be rejected at the ~ 95% level. However, from the small number of stars, it is unclear whether the two RV 
peaks correspond to approaching and receding portions of the same tidal stream, or separate substructures. 
However, it is reasonable to conclude that either (or perhaps both) of these velocities are consi stent with 
debris from the Virgo progenitor as modeled by us. We note, though, that ISesar et all (I2010al) find that 
although the two velocity groups are centered at roughly the same positions, one of th e stars in the negative- 
velocity group extends to lower right ascension. Similarly, in Figures 12 and 14 of IWatkins et all (|2009) 
the Pisces Overdensity appears to be closer at lower RA (and lower Galactic longitude) than at its peak. 
Both of these findings are consistent with our tidal stream model, wherein the debris on the lower-longitude 
side of the apoGalacticon is on the approaching portion of the stream. We conclude that it is likely that the 
Pisces Overdensity has its origin in the disruption of the Virgo progenitor. However, further observations 
and detailed modeling of both the Virgo and Pisces structures are needed to confirm this. 



6. Conclusion 

We have presented 3-D kinematics of candidate Virgo substructure members in SA 103, a 40' x 40' 
field l ocated at (a,5) 2 ooo = (178.75°, -0.59°), or (l,b) = (274. 6°, 59.1°). This work builds upon that of 



CD09L in which an orbit was presented for a single VSS RR Lyrae candidate from the SA 103 data set. We 



have obtained follow-up spectra of a total of 215 stars from the Kapteyn's Selected Areas proper motion 
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survey outlined in iCasetti-Dinescu et all (120061) . From the derived radial velocities, we isolate an excess 
of relatively high-velocity stars that are unexpected from smooth Milky Way populations. These stars, at 
100 < Vgsr < 350 km s" 1 , do not make up a kinematically cold "peak" in the RVs, but rather span nearly the 
entire high-velocity range. However, we show that the proper motions of most of the stars from these high 
velocities are consistent with their being kinematically associated. Furthermore, these high-velocity stars 
inhabit regions of the color-magnitude diagram that place them all at nearly the same distance, ~ 14 ± 3 
kpc, from the Sun. 

We identify a total of 16 candidate Virgo substructure members, from which we derive a weighted mean 
radial velocity Vgsr = 152.6 ±22.0 km s -1 and proper motion (fi a cos5,ns) = (-5.24,-0.91) ± (0.43,0.46) 
mas yr" 1 . The mean metallicity of the substructure is estimated to be [Fe/H] ~ -1.1. From the mean 
kinematics and the derived distance of 14 ± 3 kpc, we generate an orbit in a model Galactic potential. The 
Virgo substructure stars are on an eccentric orbit (e = 0.81^ 27) with peri- and apo-centers of R p = 5.6^[j g kpc 
and R a = 52^4 kpc, respectively. 

We searched for satellites or known overdensities that can plausibly be associated with the Virgo pro- 
genitor. Based on the positions and radial velocities of the derived orbit, we find that known tidal debris 
in the field of view of the ultra-faint dwarf SEGUE I match our predicted orbit in position and velocity. 
While it is intriguing to find such a close match to the predicted path and kinematics of the Virgo progenitor 
orbit, on e must explain the f act that the 24 stars making up the "300 km s -1 debris" near SEGUE I were 
found by Isimon et al. ( 201 lh to have a velocity dispersion of only 7.0 km s _1 . It is possible that these stars 
near SEGUE I, as well as many of the apparently cold stellar substructures in Virgo, are kinematically cold 
sub-clumps within the Virgo debris. Another pla usible match to our derived orbit is the globular cluster 
NGC 2419. This distant cluster was suggested by |CD09| to be associated with the Virgo substructure, pos- 
sibly as its progenitor. It is difficult to explain how the metal-poor ([Fe/H] = -2.1) NGC 2419 could have 
been the origin of the more metal-rich Virgo debris; rather, it is more likely that NGC 2419 fell into the 
Milky Way with the V irgo progenitor. We also note that there is a known overdensity of A-stars found by 
Newberg et al. d2003l) to span ~ 200 deg 2 around NGC 2419 at a mean distance of 83 kpc. This cloud of 



stars may be associated with NGC 2419, the Virgo progenitor, or both. 

Through the use of Af-body models of satellites on this rather destructive orbit, we show that the large 
spatial extent of the Virgo substructure does not arise if the progenitor of the stellar overdensity was a 
globular cluster or low-mass dwarf galaxy. Rather, the progenitor of the Virgo structure must have been 
a massive, Sagittarius-like dwarf spheroidal galaxy of mass ~ 1O 9 M . The N-body model of a Sgr-mass 
dSph on our derived orbit qualitatively reproduces the observed spatial extent, velocity spread, and estimated 
stellar mass of the Virgo overdensity. We thus conclude that the entire cloud-like Virgo substructure is likely 
the tidal debris remnant from a recently-disrupted dwarf galaxy. Our model also produces a "cloud" of 
debris near the position and distance of NGC 2419 (and the surrounding A-star overdensity), corresponding 
to the apocenter of the orbit. Another diffuse overdensity is predicted by the model near the location of the 
Pisces Overdensity, suggesting that the Pisces stars may have originated in the Virgo progenitor. 

We appreciate useful discussions with Kathryn Johnston about the expected properties of "cloud-like" 
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Table 2: Adopted Positions and Weighted Mean Kinematics of 16 Virgo Substructure Candidates 



Parameter Value 

RA (°) 178.749 

Dec (°) -0.586 

1 (°) 274.557 

b(°) 59.118 

distance (kpc) 14 ±3 
fi a cos 5 (mas yr" 1 ) -5.24 ± 0.43 

fis (mas yr" 1 ) -0.91 ±0.46 

Vgsr^s" 1 ) 152.6 ±22.0 a 

[Fe/H] -1.1 



"The Vgsr value quoted here is the weighted mean of 13 stars, after excluding the 3 likely RR Lyrae at g ~ 17. However, we note 
that if these candidate RR Lyrae stars have any large-amplitude RV variability, it seems to have no impact on our mean kinematics 
-including all 16 stars yields Vgsr = 151.8 ±20.4 km s~ l . 



